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^jqI We calculate the branching ratio of Ab A7 in the standard model using the 

^ . PQCD method. The predicted branching ratio i3(Af, A7) is about (4.3 ~ 

' 8.6) X 10~^, with reasonable parameter ranges in the heavy baryon distribution 

amplitude. This branching ratio is much smaller than those obtained in other 
hadronic model calculations. Future experimental data can provide important 
I I information on applicability of the PQCD method to heavy baryon radiative 
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I. INTRODUCTION 

Rare radiative processes involving 5 — > 57 at quark level are important for understanding the flavor changing 
structure in the standard model (SM). Exclusive radiative B decays also provide important information about the 
hadronic matrix elements where a heavy b-quark is involved. These processes being^are can also provide clues 
to models beyond the SM. There have been considerable studies on inclusive b S70, Q, and exclusive mesonic 
B K*^^ both experimentally and theoretically within and beyond the SM 1^,1^ IS- Theoretical predictions for 
inclusive decays agree with data very well in the SM. Calculations for exclusive processes are in general consistent 
with data although there are unavoidable uncertainties due to our lack of good understanding of QCD at low energies. 
Nevertheless methods have been developed to calculate hadronic matrix elements in recent years 0, |^ . With more 
data becoming available, new b-decay processes can be studied. These processes can be new tests for different methods 
in calculating hadronic matrix elements and new physics beyond the SM. In this work we study Afc A7. In this 
decay more experimental information about the heavy b quark inside the hadron which is not available in inclusive 
and mesonic b-hadron decays, such as spin polarization during hadronization, and the handedness of the couplings 
at the quark level, can be extracted(9t lldLllllll2Lll3j . Therefore the baryonic b-hadron radiative decay can provide a 
new test for theoretical methods for b-quark hadronization. 

There are some studies in the literature on A7 0,0, 0,0,0 decay ranging from phenomenological models 

to QCD sum rule approaches. Our study will be based on the PQCD method [l^ lisl IT6|. This method has been 
shown to give consistent results for two body mesonic B decaysQ. We expect a PQCD calculation for Ah — > A7 will 
also give a reasonable estimate since the energy-exchange carried by gluons in the matrix element calculations is large. 
Result obtained in this way can serve as a good reference for discussing the relevant hadronic matrix elements. 

For SM, the effective Hamiltonian responsible for b s"/ comes from the electromagnetic penguin diagram and is 
given by [131 : 

Heff - i^Vt,V*£^C',"{fi)m,-s<j^^{l+j,)bF'^'^, (1) 

where C^^^ = mi,) = —0.31. In our numerical calculations, the running of Cy''^'^ will also be taken into account. 

It has been shown that there may be resonant (long distance) J/'0(V'') contributions . If these contributions are 
included, one should add a term (3Ci(^) + C2i^i)){^/al^) Ej^v.'A' ^j(0)%'^r(j ^ l+l-)Mj/{q^ - Mf + iMjTf*) to 

the Wilson coefficient C^-^-^ . Since for 6 ^ 57 process, — 0, there are double suppressions for the long distance 
resonant contributions with one of them coming from the Breit-Wigner factor ~ TjJMi and another coming from 
the extrapolation ofujj{Mf ) = 1 to Wj(0) with uij{0) < 0.13 (and could be smaller) [isj . we will neglect the resonant 
contribution for radiative decays in our later discussions. 

At the hadron level, the decay amplitude for Af, — > A7 is obtained by inserting the effective Hamiltonian between 
the initial and final hadron states, 

M(Ab-. A7) = (A7|i/e//|Ab). (2) 
There are two form factors for A^ — > A7 from the above which we write as 

M,, = {A{p')\C^"{fi,0)sa^,q-'{l + 75)6|Ab(p)) - M^iFLa^^q" [1 - 75) + Fna^^q^il + i^))Kb{p). (3) 
We obtain 

r(A, ^ A7) = ^^'^*^'''^/3^-J^^^''''"^' « - mlf{\F,\- + \Fnn (4) 

Emission of a photon from the tree operators O1.2 can also contribute to A;, A7. Although the Wilson coefficients 
of these operators are larger than those of the penguin operators, there is a large suppression coming from the CKM 
factor |Kif)l^*5/Vt{,Vt* |. The overall contributions from bremsstrahlung of a photon off the operator O1.2 is therefore 
suppressed. We will neglect their contribution in rest of discussions. 

II. PQCD CALCULATION OF THE HADRONIC MATRIX ELEMENTS 

We now describe our calculations for the hadronic matrix elements defined above using the PQCD method developed 
in Ref. [T3 . ll5lll^ . We define, in the rest frame of A^, p, p' to be the Af,, A momenta, ki{i = 1, 2, 3) to be the valence 
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quark momenta inside Af,, and k[ to be the valence quark momenta inside A. We parameterize the light cone momenta 
with all light quark and baryon masses neglected as 

p^ip+,p-,OT) - ^(1,1, Ot), p' - (y+,0,OT) 

fci = (p+,2:ip",kiT), k2 ^ iO,X2P^ ,'k2T), /ca = (0,X3p~,k3T) 

k[^ix[p'+,0,k[j,), fc^ = (xip'+,0,k^j,), fc^ = (x^p'+,0,k^j,) (5) 

where Xi and x'^ are the fractions of the longitudinal momenta of the valence quarks with xi + X2 + x^ = 1 and 
x'l + x'2 + x'^ = I. kiT and k^-j^ are the transverse momenta of the valence quarks inside A;, and A, respectively. 

As a self-consistent check, one should make sure that the expected relation — fcj ~ 0{-^QCDmt) holds, since the 
light quarks in the heavy baryon should have momenta of order Aqcd- Naively, the above gives a value of order 
(1 — ^2)™^^ which does not have the explicit form as expected. To understand this, one needs to combine the form 
of the heavy baryon wavefunction which determines how quark momenta are distributed inside the baryon. We have 
checked this using the wavefunction given later, and obtained the ratio of average values < X2.3 > / < xi >~ m/mA,,, 
where m is of order Aqcd- With the constraint xi + X2 + X3 = 1, the desired order for p^ — is then obtained. 

One can write = pp~^ with p = ^f- = ^— fe— The ranges for and p are given by < < (Ma^ — m\)'^ 

and 2m!y/Mf^^ < P < (-^-^A^ + "^DZ-^A^ off-shell photon is allowed. In our case of A^ A7, = 0. Here we have 
kept A mass in the expressions for the purpose in tracing the ranges of the kinematic variables. In the approximation 
we are using, it should be set to zero as mentioned above. 

In the PQCD picture, hadrons are formed from quarks with appropriate wave functions describing the momenta 
distribution of quarks inside the hadron. The Af, wave function is usually defined through the quantity [l9l |20( . 

= ^^[{i> + Mf,,)^^C]p-\K,{p)]a.^{h,v), (6) 

where /a^ is a normalization constant, Af,(p) is the A;, spinor, and ^(ki,fi) is the wave function. Here we have used 
the heavy quark symmetry which should be applicable in the present case, following Refs.^3|23|, to reduce the form 
factors to the above simplified form. In general there are more components in the wavefunction if all quarks are light. 
For the light baryon A the leading-twist wave function of A is defined byp^: 

2 _ 

(yA)„^^(fc» = J^^:!^01e^K^'e^'^^O\TK{w^)uliw2)d^O)m^^^^ 



-{(j/C)^^[75A(p')]a<i>^(fc» + (;/75C)^^[A(p')]a$^(fc»} 



~^^{a,,p'-C)p,[ri5A{p')U^''{kl u), (7) 

where /a and /a are normalization constants, and A(p') is the A spinor. 

Including the Sudakov factor with infrared cut-offs w(w'), and running the wavefunction from v down to w(ijj'), 
then we obtain[l3|: 

*(a;i,6i,p, = exp[-^s(w,a;;p") - 3 / —-iq{as{jl))]^{xi) 
1=2 ^ 

<i>^x',,b',,p',,,)^exp[~J2s{uj\x[p+)-3 r ^j^iaMWix',), 



(8) 

|bp-b 

and 6^ ' = |bi ' - ^ |. Here b and b' are the conjugate variables to ky and kj. defined in Appendix B. 



where j = V, A, T,uj = mm(l/6i, I/62, l/h), and cu' = min{l/b[, l/b'2, V^s)- = - bp|, 6^ = |bp - b^'^|. 
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The explicit expressions for the Sudakov factors are given in Ref. 14'| with 

s{u;,Q) = p ^[lnfi)A[asip)]+B[asip)]] 

J UJ P P 

^ as r67 TT^ 10 8^ , ,e^^ ,as 
A = C^.- + [y-y-27"/ + 3/3oM-)](-)^ 

3 TT ^ 2 ' 

Po ^ (9) 

where 7^ is the Euler constant, is the flavor number, and jq is the anomalous dimension. For Af, baryon decays, 
the typical energy scale is above the charm mass. We will take n/ equal to 4 in our calculations. 
The hadronic matrix elements can be written as: 



M,,^ = J[Dx] j[Db]{YA)^,p,,,{xi,bi,p',,.) 

HtJ "'^'^ , a; ■ , 6^ , 6 ■ , Ma, , i^) (Fa J (x, , 6j , p, J/) , (10) 

where the measures of the momentum fractions |l4j are give by 

3 3 

[Dx] — [dx\[dx'], [dx\ = dxidx2dx3S{l — xi), [dx'] = dxidx2dx'^S{l — ^^xj). (11) 

1=1 1=1 

The measures of the transverse extents [Db] are defined in Appendix A. 

The hard scattering amplitude H"^ {x, x' ,b,b' , Ma^,^) is obtained by first evaluating the amplitude 

^ "'''''(xi, x^, k-r, k^, MaJ for the 'i'th diagram in Fig. for a corresponding Wilson coefficient Cf''^'^ which 

is displayed in Appendix B. One then carries out a Fourier transformation on k^ and k(p to b and &' space to obtain 

H;'^^''^'^'"'^^{x, x\ 5, 6', A/aJ. The procedure of carrying out this transformation is described at the end of Appendix 
B.' 

Collecting all contributions in Fig^and multiplying the corresponding Wilson coefficients, one then obtains a hard 
scattering amplitude H^J^ ^ "^'''(x, x' ,b, b\ MaJ — J2i ^i'^'^ ^ ^ "^^{x, x' , b, b' , A/aJ- Here we have labelled the 
hard scale as t which is taken to be the larger of the two variables ii_2 associated with the virtual gluon momentum 
in Fig. n i.e. t = max.{t\^t2). The expressions for ti^2 are listed in Appendix C. 

Finally a RG running is applied to the hard scattering amplitude to match the scale v in the wave functions and 
we obtain (1^ 

Hf^P'^''^P'<{x,x'Ab',M^,,v) = 



exp[-6 / " ^7,(«.(7^))] X Hlf^'^-^^x, x\ b, b', MaJ (12) 

The form factors are obtained by grouping relevant terms according to the definition in eq. Using ea. pUf) we 
obtain a generic expression for the form factors corresponding to each diagram as 



^fih, I [DA f met" if 

j=V,A.T 

^ ^,{x)^\{x')<iM-S'W^^' 

3 3 /-t Jtt l-t 



S = ;^s(a;,Xfcp-) + ;^s(^',4p'+) + 3 / ^7,(«.(m)) + 3 / ^7,(«.(P)) (13) 



k=2 k=l 



where represents the form factors contributed by the "i" the diagram in which operators with the Wilson coefficients 
Cf" are inserted, in our case Cf'^^ = C^" . The superscript j labels V, A, and T related to the spin structure of the 
valence quarks in the A baryon with = = /a. The explicit expressions of are presented in Appendix D. The 
functions Hp are given in Appendix E. The total form factors are obtained by summing over contributions from all 
diagrams. 
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FIG. 1: The lowest order diagrams for the A;, — + A7 decay. The sohd hnes, double lines, wavy lines and the black blub vertex 
denote the light quarks, b quark, gluon and the electromagnetic penguin vertex, respectively. Diagrams with triple-gluon vertex 
do not contribute since their color factors are all zero in the present case. 



III. NUMERICAL RESULTS 



We are now ready to evaluate the form factors numerically. For concreteness, we adopt the model proposed in 
Ref.^3 fo'^ the baryon distribution amplitude \E', 

M? 

.3) = iVx,x..3exp[-^ -Wk-Wtj- ^^'^ 
The normalization constant N is obtained by the condition: 



J[dx]^{xi,X2,X3) ^1. (15) 



The decay constant /a^ is determined by fitting B{Ab — > AJV) whose central value is 5% measured by DELPHI |21| 
using the same PQCD method. When fitting the data we truncate the double log Sudakov factor in such a way 
that the factor exp(-s) is smaller than 1 following the prescription in Ref.'22'|. Our numbers for /a^ are different 
from those obtained in Ref. 15] where a B{Ai, AJV) was taken to be 2%. We also have chosen cut-offs as 
uj = 1.14mzn(l/&i, 1/62, 1/&3) and uj' = 1.14mm(l/6']^, I/62, 1/63). The factor 1.14 is adopted because this cut-off 
choice can result in form factors which vary smoothly with square of momentum transfer in fitting A^ — > AJV process 
and it reflects the resummation of next-to-leading double log in higher order radiative corrections 23J . Also the f3 and 
ruq in the heavy baryon wavefunction distribution need to be fixed. In Ref. [l3 . Iisl IT^ . /3 = 1 GeV and = 0.3 GeV 
were used to estimate Ab — > AJD, Aj, plD and also Af, — > AJ/ip decay rates. f3 should not be too much smaller 
than 1 GeV if the form factors are dominated by perturbative contributions. Therefore we will let both P and rrig 
vary within ranges as 0.6 ^ IGeV and 0.2 ~ O.SGeV. The results for /a,^ are shown in TableQlfor different parameter 
choices respectively. 

The A baryon distribution amplitudes have been studied using QCD sum rules. In this work, we adopt the model 
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TABLE L Decay constant /a^ for different choices of /3 and m,, respectively. 



/A.(GeV) 


13 = 0.6GeV 


P = 0.7GeV 


P = O.SGeV 


P = 0.9GeV 


P = IGeV 


ruq = 0.2GeV 


0.691 X lO-" 


0.841 X 10"^ 


1.02 X 10"^ 


1.21 X lO"'' 


1.43 X 10"'' 


niq = O.SGeV 


1.27 X 10~^ 


1.45 X 10"'' 


1.65 X 10"^ 


1.88 X 10"^ 


2.12 X 10"'' 



proposed in Ref. [2J| , 

(l)^{xi,X2,X3) = A2(f>asixi,X2,X3)[0.18{xl - xj) + 0.10(a;2 - 2:3)], 

(l3^{xi,X2,X3) = -42(/)„^(a;i,a;2,a;3)[0.26(a;^ + +0.34a;2 - 0.56a;2a;3 - 0.24x1 (2:2 + 2^3)], 
(j)'^{xi,X2,X3) = 42(/ias(a;i,X2,a;3)[l. 2(2:2 - X3) - 1. 4(2:2 - 2:3)], 

0Qs(xi, 2:2, 2:3) = 1202:1x2x3. (16) 

The asymmetric distribution in the momentum fractions of the three quarks impHes SU (3) symmetry breaking. 
The constants /a and /J are fixed to bep^l 

/a = 0.63 X 10"2GeV^ /J 0.063 x lO^^GeV^. (17) 

Finally to obtain the branching ratio for A;, — > A7, for definitiveness we fix rest of the parameters as following. The 
parameter Aqqd which enters in the strong coupling constant and various Wilson coefficients, the b quark mass and 
the CKM mixing parameters are set to be: Aqcd at 0.2 GeV, mt, — 4.8 GeV, and the CKM mixing parameters are 
set to their central values'^: S12 = 0.2243, S23 = 0.00413, S13 = 0.0037 and ^13 = 1.05. 

Our explicit calculations show that Fl = and a non-zero value for Fjf as expected since light quark and light 
baryon masses have been neglected. The contributions from each diagrams for Fr are shown in Appendix E. The 
resulting branching ratio is shown in Table |nj We see that the branching ratio for A;, A7 is in the range of 
(4.3 ~ 6.8) x 10-8. 



TABLE IL Branching ratio(BR) of Af, A7 for different choices of P and with Sudakov truncation. 



BR(xlO*') 


P = 0.6GeV 


P = 0.7GeV 


P = O.SGeV 


P = 0.9GeV 


P = IGeV 


m, = 0.2GeV 


6.76 


6.26 


6.19 


4.90 


4.67 


rUg = O.SGeV 


6.42 


5.75 


5.61 


4.44 


4.32 



IV. DISCUSSIONS AND CONCLUSIONS 



In this work we have used the perturbative QCD approach to evaluate the branching ratio for radiative decay 
Afe — > A7. This process occurs via penguin diagrams. Our results are shown in Table ITU The branching ratio obtained 
is much smaller than results obtained, shown in Table UTTl using other methods. 



TABLE III: Decay branching ratios (B) of A;, —* A7 based on the form factors from the QCD sum rule approach, the covariant 
oscillator quark model, HQET and MIT bag model, respectively 



Model 


pole modelflO] 


QCD sum rulefll] 


covariant oscillator quark modelfl2] 


HQETflS] 


bag modelfl^ 


B 


(0.10 ~ 0.45) X 10"" 


(3.7 ±0.5) X 10"" 


0.23 X 10"" 


(1.2 ~ 1.9) X 10"" 


0.6 X 10"" 



There are uncertainties in PQCD predictions due to unknown parameters in wavefunctions. We have tried to under- 
stand such uncertainties by varying several relevant parameters. Within reasonable ranges of the parameters it is not 
possible to obtain a branching ratio larger than 10~^. We have considered another possible uncertainty in the method 
used here. This is the choice of the infrared cut-offs lu{ll!') in the Sudakov suppression factor which damps the perturba- 
tive contributions. In our calculations the cut-offs are set to the conventional values with uj = 1.14rnm(l/6i, I/62, 1/&3) 
and uj' = 1.14min(l/6i, I/62, 1/^3) discussed in the text. The factor 1.14 is adopted because this choice for cut-offs' 
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can result in form factors which vary smoothly with square of momentum transfer in fitting Aj, AJV process and 
it reflects the resummation of next-to-leading double log in higher order radiative corrections [23 ■ We have checked 
with slightly different cut-offs and find impossible to obtain branching ratio to be as large as what listed in Table ITTll 
The prescription of truncating the factor exp(-s) to be smaller than 1 described in Rcf.'22'| may also be a source 
for uncertainties. We therefore have evaluated the branching ratio without this truncation. The results are shown in 
Table ITVI We see that the results are similar to those obtained in Table ITTI 



TABLE IV: Branching ratio(BR) of At A7 for different choices of (3 and without Sudakov truncation. 



BR(xlO'*) 


P = 0.6GeV 


P = 0.7GeV 


(3 = O.SGeV 


(3 = 0.9GeV 


/3 = IGeV 


m, = 0.2GeV 


8.60 


7.22 


5.91 


4.92 


4.60 


= 0.3GeV 


5.96 


5.73 


5.70 


4.67 


4.30 



We therefore conclude that within the PQCD framework, the branching ratio for Af, A7 is much smaller than 
other model calculations. This is somewhat surprising since PQCD calculation for the branching ratio of B ^ K^*'>j 
obtains a value of order consistent to other model calculations and also agrees with experimental value of about 
4 X 10~^[^. There is a huge suppression for Af, — > A7. At this moment there is no data available for A;, A7 yet. 
One has to wait for future experimental data to tell us more. If a branching ratio above 10"'' is measured at some 
future facilities, such as LHCb, the PQCD method used here will certainly need to be modified. 

On the theoretical side, one expects the branching ratio for Af, — > A7 to be smaller than that of B ^ K^*'>j due to 
several suppression factors such as an additional and a large momentum squared suppression factor as one more 
hard gluon is exchanged between quarks. There is also an additional Sudakov suppression factor due to an additional 
spectator quark involved in the process as can be seen from ea. (|13|l . 

One might question the applicability of PQCD method for the process under consideration. One notes that in 
the PQCD approach, both gluons are hard ones which excludes the possibility of including contributions where two 
spectator quarks (not involved in the weak interaction vertex) form a collective object first due to soft gloun exchanges, 
i.e. the diquark, and then this object interacts with the other quark by exchanging a hard gluon. If this contribution 
turns out to be the dominant one, the branching ratio may be substantially larger. At present there is no solid 
theoretical method to treat this effect yet, we do not have a definitive answer about this. We, however, note that 
estimate for Af, AJ/ ip using the same method gives a reasonable range compared with data[l2 . This can be taken 
as a support for the applicability of the method to Af, decays. Our result for i?(Af, — > A7) represents a reasonable 
estimate. The branching ratio for Af, A7 is in the range of (4.3 ~ 8.6) x 10^^ which is smaller than predictions 
using other methods listed in Table Hill We have to wait for future experiments to provide more information. 
Acknowledgements : 
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Appendix A: the b measures 

The ordinary b measure is defined as 

m = ^ (18) 

The explicit forms of [DhY for each diagram i in Fig. 1 are given by 



[Dh] 


a) ^ 


[dbi][db3][db'i][dby. 


[Dh]^''^ 


= [dhi][dh2][dh[][dh'2 


[Dh] 




[dbi][db3][rfb'i][dm. 


[Dh]^''-^ 


= [dhi][dh2][dh[][dh'2 


[Dh] 


= 


[db2][rfby[rfby. 


[Z3b](-^' 


= [db3][dby[dby. 


[Dh] 


g) = 


[db2][db3][db^5]. 


[DbjC*) 


= [db2][db3][dm 


[Dh] 


*) = 


[dhi][dh2][dh[][dh'^], 


[Dh]'-J^ 


= [dh,][dh3][dh[][dh',] 


[Dh] 


fe) ^ 


[dbi][db2][db;][dm. 


[Db](') 


= [dhi][dhs][dh[][dh's 


[Dh] 


m) _ 


- [db2][db3][dby. 


[Db]("' 


= [dh2][dhs][dh's]. 



Appendix B: Hard scattering amplitudes H^'^ ^ ^ "''^(x^, x^, ky, k^, MaJ 
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Expressions of amplitude H^'^ ^ "'^'^ {xi, x'^,\^T^^T^ ^A^) for each diagram in Fig. 1. In the following Oj^ comes 
from the 7-matrix in the effective Hamiltonian, O^^^ = a^i^q'^ R 



For the hard amplitude of Fig. 1(a): 

H-/0'^'-P^{x,,x[,\^T,^^T.M^,)^ 

\e £ {l-')c'c{J--' 1 )b'b{l )a'a\gs 



ip' -k[- k3)^{p' -p + ki)^p-p' + k[ - kiYik^ - k'^Y 

"^^^ [Aa + {Kt + k3T)2] [Ba + k^j,] [Ca + (kiT - KtY] [Da + (kgT - KtY] 



with 



Aa = X3(l - x[)pMl,Ba - (1 - xi)pMl,Ca = (1 " a^i)(l - x[)pMl,Da = xsx'^pMl (21) 
and the color factor 

For the hard amplitude of Fig. 1(b): 

^b,a'/5Va/37(^^^^^^kT,k^,MAj = 

rabc^a'b'c'(rJ.^rJ.j. /rj.-,^ frj^j. 14 {Y)l3'p[lpiJ>' ~^'l~h)l\]l'l[l^{f -i>+h)0,,]a'a 
[e e (1 l-')c'rM )b'b[J--' )a'a\gs 



ip' - k[ - k2)Hp' -p + kiY{p -p'-ki + k[Y{k'^ - k2\ 



2 



with 



[Ab + (k'lj, + k2T)2] [Bb + kir] [Cb + (kiT - Kt?] [Db + (k2T - k^T)'] 



Ab ^X2{1~ x[)pMl,Bb = (1 - xi)pMl.Cb = (1 - xi)(l - x[)pMl.Db = X2x'^pMl^ (24) 



Inspection of the above calculations, one notices that one can easily obtain i?^'" ^ "'^'''(a;^, a;^, k^, k^, Ma,^) from 
p 7 "/37(2;^^ 2;^^ if^, MaJ and vicc versa by simply exchanging the momentum indices 2 and 3 for k and k', 
and exchanging the positions of the Dirac indices j'j and Due to these properties, the contributions to the 

form factors from the above two diagrams are the same. This fact can be easily understood by noticing the following 
properties of the quantities related to the distribution amplitudes: i) The distribution amplitudes '^(xi,X2,X3), and 
(f>^{xi, X2,X3) are symmetric in exchanging X2 and xs, while (f>^''^{xi, X2,X3) are anti-symmetric in exchanging X2 and 
xs, as can be seen from eas. (|14|l and (|16|l . And ii) When exchanging the Dirac indices /3 and 7, the expressions for 
{YA^)ai3-y{kii') in eq.lO, and terms proportional to (j)^ for {Y\)aj3j{k'i, f) in eq.jTI) will have a sign change, while terms 
proportional to <j)^''^ remain the same. Since going from the contribution of diagram (a) to diagram (b) involves both 
actions: exchanging the momentum indices 2 and 3, and the Dirac indices f3 and 7, this results in no sign changes for 
all the terms involved. After integrating out x{x')2,3 and 6(6')2,3 to obtain the final form factors using ea. (|13|l . one 
then obtains the same results for both diagrams (a) and (b). 

Similar situation happens for the following pairs of diagrams: (c) and (d), (e) and (f), (g) and (h), (i) and (j), (k) 
and (1), and, (m) and (n). In the following we will only display the results for diagrams (a), (c), (e), (g), (i), (k) and 
(m). The expressions for diagrams (b), (d), (f), (h), (j), (1) and (n) can be obtained by exchanging x(x')2 and a:(a;')3, 
and also 7'7 and /3'/9. 

For the hard amplitude of Fig. 1 (c) : 

^c,o';3Va/37(2,^^^^^kT,k^,AfAj = 

rabcWb'c'irj.j-, tT^trpn i 4 ilp)^ - h ~ ^3h^]f3' pilxif - ^ + h)0 ^,U' a 

[e £ {l-'jc'cK-l- J--')b'b{^ ja'aWs 



{p~ki~ k'^Yip' ^p + kiY{p -P' + k'l - kiYih - k'^Y 



with 



"^^^ [A, + (kiT + k'sTY] [Be + kl^] [Cc + (k3T - k'^^Y] [Dc + (kiT - K^Y] 

A, = x'^{l~ xi)pMl,B, - (1 - x,)pMl,Cc = xsx'^pMl^D, = (1 - xi)(l - x[)pMl (26) 



For the hard amphtude of Fig. 1(e): 

{e^'"e'''''''{r)e,4rT^)b'biT')a'a]9. 
= CN9t 



{p' -k'^- h)^p' -p+ fci)2(fc^ - fc2)2(A;3 - 

(7p)7'7[7''(/ -1^2- h)iHf -i> + h)o,U'M0'p 



with 



= a;3(l - x[)pMl^,B^ = (1 - xi)pMl^,Ce = X2x'^pMl^,D^ = x^x^Ml^ 
For the hard amphtude of Fig. 1(g): 

H9/0'i'-Pi{xu X\, kT, k^, MaJ = 

^^ahc^a'h'c' ,rj.i^ (rj^ys frptj^j. i„4 (7p)7'7[7^(/ " 1^2 " h)0 t^{i> ' h " 1^2 + '7^b)7^]a'a(7A)/3'/3 
[£ (i Jb-bU [(p_^^^_fc,)2_^2](p,_fc,_fc3)2(fc,_fc^)2(fc._fc3)2 

(7p)7'7[7^(/ ^ 1^2 - h)0^^{i> -h-^2+ r^b)7^]a'a(7A)/3'/3 

[Ag + (k'^r + k3T)2][5g + (k'^r + k3T)2][Cg + (k2T - KT^iDg + {k^T " K^f] 



CN9t 



with 



Ag = (x^(l - X3)p + X3)Ml^,Bg =X3{1- x'2)pMl^,Cg = X2x'2pMl^,Dg = xax'spMl^ 
For the hard amphtude of Fig.l(i): 

^i,aW«/37(a;.,a;^,k^,k^,MAj = 

^'(T*T^)c'c(r)b'b(T^)a'a]fff 



-ft[- hhxh'^jOAi^ + mh^'ahn^'ff 

[{p -p' + fci)2 - ml]{p' -k[- k2f{p -p'-ki+ k[f{k'^ - k2f 

n 4 [7p(/ - ^1 - h)l>\^"^{0^{i> -i>'- + v^h)l^V'o.{Y)?'P 

""^^ \Ai + (k2T + Kt?\ \Bi + (k'lr + k2T)2] [C, + (kiT - k;r)2] [I), + (k2T " ^2T?\ 



with 



= (1 - x'^)pMl^,Bi = X2{1 - x[)pMl^,Ci = (1 - xi){l - x[)pMl^,Di = X2x'^pMl^ 
For the hard amphtude of Fig.l(k): 

Hl''^'0'^''^P'^{xux[,kT,^^,Mt,,) = 

^^abc^a'b'c' (rp^rr.j. frpj. frJ.^^ .4 [7p " ^1 " ^2 )7a] 7'7 [^^ " 7^' + l^'l + "^6)7l «'« (7^ )/3'/3 

[e e (J l-')c'c{-l--')b'b{-L )a'aj5s 



[(p - p' + fci)2 - m2](p - fci - k'^Y{p -p'-ki+ k[y{k'^ - fc2)2 



[A, + kf^][S, + (kiT + k;T)'][t^ft + (kiT - k'ir)2][l?, + (k2T - k^r)"] 



with 



Ak = {l- x[)pMl^,Bk = 4(1 - xi)pMl^,Ck = (1 - a;i)(l - x[)pMl^,Dk = xz^P^i 
For the hard amphtude of Fig.l(m): 

^m,a'/3'7'a/37(^.^ ^/^ ^t, k^, Ma J = 

!.abc.a'b'c'(rpj^ /^iN (TiT3\ 1^4 (7A )-, '-.J*^^// (/^ - / + l^'l + >l>b)-pifi-h - f'j + 'U..'nA^; 

[S (1 Jb'6U VaJ,9. [(p_p,^^)2_^2][(p_^^_^,)2_„2](^,_fc^)2(fc, _^3)2 

4 (7A)7'7p/. (:?^ - / + K + rnfc)7p(?^ - - ^.^ + mbh%,c,{-fP)[rp 

""^^ [Am + k'^ABm + (k2T + KTr][Cm + {^2T " Kr^lDm + {^3T " k^r)"] 
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with 

Am = {l- x[)pMl^,Bm = (4(1 - X2)p + X2)Ml^,Cm = aJs^pM^^, = xax'spMl^ (36) 

The expressions for the hard scattering amphtude in b and b' space are obtained by making a Fourier transformation 

on and k'j, space. In the following wc given one example for Fig. 1(a) as an illustration. We note that the kx and 
k'j, dependencies are all in the denominators in the above expressions, one then just needs to consider that part of the 
fourier transformation. For Fig. 1(a), it is given by 

n'-''\xi,x',,kT,k'T,Mj^,) = -, ^ ; -, . (37) 

^ ' ^ [A„ + (k;r + k3T)2][Ba + k?^][a + (kiT-k;y)2][£'a + (k3T-k^T)'] ^ ' 

The fourier transformed expression is then given by 

n^''\xi,x'i,bi,b'„MA,) = j e-^(''^-^^+'''i--^'i+''^-''^+''3T-^)n(a)(a;.,a;^,kT,k^,MAj(i2kird2k;^^^ (38) 
Defining k^r = + ksT, ^bt = kiT, kcT = ^it — k^j., and kjjT = ksT — kg^, we rewrite the transformation as 

n^-\x„xi,b,,b'„MA,) 

g-i[kAT-(b3+b^)+kBT-(bi+b'i-b3-b^)+kcT-(-b'i+b3+b^)+kDT-(-b^)]X 



-I 



k^T c?^ ksT c?^ kcT rf^ k/jT 



(klr + Aa)(kj,T + Ba)(k'cT + CaXKx + Da) 

(27r)4/fo(v^|&3 + b's\)Ko{^/K\bi + b[ - 63 - b's\)Ko{VC^\b[ - 63 - b's\)Ko{^/D~a%\), (39) 



In the above we have used 

gjkb 



d^fcp-p-^ = 2TrKo{y/A\h\), A>0. (40) 
One obtains the expression for H^'" ^ i" °'f^'<{xi, x[, b, b', A/aJ as 

H-'^'P'y"f^^Xux'i,b,b',MA,) = CNgtilp)y-y[Y{f -^[- hh^h'0bx{f - ^+ hWa'a^Hxi,4,bi,b'„MA,). 

(41) 

In carrying out the fourier transformations for other diagrams, two other forms of functions will be encountered. 
We list them in the following 



/ 



(fc2 + A)(p + B) Jo y/z; 

where ^4 > and B, C arbitrary. Ko and Ki are the modified Bessel functions of the second kind. And 

Zi = Az + B{\- z), 

Z2 = A{1 - Z2) + {B(\ - zi) + Czi], 

^i(l - Z\) 

X2 = (bi - z^h^f + ^i^i^b^ (43) 

Appendix C: The maximum of 2 

The hard scales, the maximal of t\ and t\ for diagrams (a), (c), (e), (g), (i), (k), and (m) in Fig. 1. Exchanging 
b{b')2 and b{b')z, one obtains the expressions for diagrams (b), (d), (f), (h), (j), (1) and (n). The expressions of Ci,Di 
are collected in Appendix B. 
Appendix D: Expressions of fl^ 

The expression of fi* for diagrams (a), (c), (e), (g), (i), (k), and (m) in Fig. 1. Exchanging b{b')2 and 6(&')3; 
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(a) 
(c) 
(e) 
(g) 
(i) 
(k) 
(m) 



max{y|Cc|, 1^ 
max{\/|Cel, 7E^i^>'^'} 

1 



max{^|Ci|, 
max{-y|Cfe| 
max{^/|C7^ 



|bi-b2-b^ 



S 

max{\/|Da|, -pgrj ,tj,^i;'} 

max{y|_Dc|, i^.tJ.cj'} 

max{v'|_De|, 

ma.x{^\Dg\, j^,uj,Lo'} 



"3 

max{^jA|, 

UJ,U}'} 



max{\/|_Dfc|, 
max{^/\D^\, 



1 

lb2l ' 
1 

Ibal ' 



obtains the expressions for diagrams (b), (d), (f), (h), (j), (1) and (n) 



with 



n^-^ = (2^)4Xo(VX|b3 + h'^\)Ko{y%\h, + b'l + b3 + h'^\)Ko{Vc'c\h[\)Ko{y^c\h3\) 

Jo ^lU ^ ^l) V ^^2 



n(3) = 167^5 / dz 



|b3 + b^|/fi(yzf|b3+b: 



^Xo(v^|b2|)i^o(v^|b:,|) 



f7« = (27r)4Xo(VA;|bi + b'l - b2 - m)i^o(/B;|b2 + b^Di^oCVQIbiDifol/D^mi) 
n^''^ = (27r)4i^o(V^|bi + b'l - b2 - b'2|)i^o(V^|b2 + m)i^o(VC^|bi - b2 - h'^\)Ko{^/D^\h2\) 

1 



dzidz2 



^i(l-zi) v^Z^ 



ifi(v/X™Z^)i^o(v^|b3| 



X| = (b^ + Zib, 



)' + ^i^^^bl,Z|=A,(l-Z2) 

Z2 



Z2 



Zl(l - Zl) 



Z{=AgZ + Bg{l-Z) 

2 , - 2l) 



X™ = (b'2 + Zib2 



Z2 



^2 ' ^2" 



A,n{l-Z2) + 



Z2 



Zl(l - Zi) 



[Be(l-Zi)+CeZl] 



[B,„(l-zi) + C'„zi] 



(44) 



(45) 



Appendix E: Expressions for Hp 



In this appendix we list Hp corresponding to the form factors defined in ea. (|13|l . We use for each diagram. The 
expressions for diagrams (a), (e), (g), (i), (k), and (m) in Fig. 1, whenever non-zero, are hsted in the following. The 
expressions for diagrams (b), (f), (h), (j), (1) and (n) can be obtained by exchanging x{x')2 and x{x')3 and changing 
the signs for expressions F^'^ ■ Diagrams (c) and (d) have no contributions to A;, A7. Fl is equal to zero in our 
approximation. 

For the hard amplitudes of Fig. 1(a): 



F^ = 8x3pH4 



(46) 



The relation between the tilde form factors listed above and the form factors in eq.|(2l is as the following, taking 
F^ as an example, F^ = f^/i/^^ J[Dx] J[DbYCf^ {t')^'A,{x)^i{x')exp[-S']F^n\ For this example j = A, and 
/a = /a- For Fig. 1(a), 'i' takes the value 'a'. Similar for other form factors and diagrams. 
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The other non-zero contributions are 

Fig.l{e) : = ^ -U^yx^, 

Fzg.lig) : = ^Mlinibp + MA,(-2a;3(-l + p) + (1 + + P))p)) 

= 4Ml^(-l + p){mt + Ma,(-1 + 

Fzg.l(fc) : Fi^8Mlpx'^{mb + MA^i-l + p)), 

Fig.l{m) : = -4M|^(-m2(-2 + p) + MA,mb(a;2 - + (-1 + x[ + x'^)p) 

+Ml(-2 + (2 - x[ + x',)p - ^ ^^(^ ^ ^ ^, 

-4M|^(mgp + MA,mf,(-2 + 0:2(1 + p) + (-1 + + 4)p) 
+M2^(2 - (2 + x'l + 4)p + + + (1 + ^'^)^))) (47) 
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